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Abstract 

Carrying digital information in traditional copper wires is becoming a major issue in electronic 
circuits. Optical connections such as fiber optics offers unprecedented transfer capacity, but the 
mismatch between the optical wavelength and the transistors size drastically reduces the coupling 
efficiency. By merging the abilities of photonics and electronics, surface plasmon photonics, or 
'plasmonics' exhibits strong potential. Here, we propose an original approach to fully understand 
the nature of surface electrons in plasmonic systems, by experimentally demonstrating that surface 
plasmons can be modeled as a phase of surface waves. First and second order phase transitions, 
associated with percolation transitions, have been experimentally observed in the building process 
of surface plasmons in lattice of subwavelength apertures. Percolation theory provides a unified 
framework for surface plasmons description. 
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Surface plasmons (SPs) plasmonics is of interest for a wide range of domains, from physics, 
chemistry, to biological sciences . SPs are associated with electromagnetic waves trav- 
eling along a metal/dielectric interface with collective motions of electrons. The tight lo- 
calization of SPs at the surface of the plasmonic circuit provides the possibility of focusing 
and waveguiding light in subwavelength systems, which can closely interact with electronic 
devices. In the recent years, a renewed interest in SPs has arisen from the demonstration of a 
strong and unexpected enhancement of light transmission through arrays of sub-wavelength 
holes [3|. Enhancement of several orders of magnitude has been reported with respect to 
standard aperture theory 0, 0, |f|. The transmission can even exceed the surface ratio occu- 



pied by the holes, implying that light is focused by the structure o 



the arrays through the 
10] : near field microscopy, 



holes. These results are stimulating in numerous fields 
high density storage, detection of chemical and biological molecules, light emission, optical 
nonlinear devices. 

Since this discovery, many studies have been performed to fully understand and control 
SPs, involving the influence of the shape, geometry lattice and polarization in one- and 



two-dimensional subwavelength structures U|, [12|, ll3|, |14|, ll5[. However, the fundamental 



question of the establishment of the SPs themselves has never been fully addressed yet. 
Here, we experimentally demonstrate that one- and two-dimensional SPs can be generated 
through a phase transition like process, separating phases of surface waves, i.e. specific 
organizations of the electromagnetic waves at the surface of the perforated metal plate. The 
analysis of the phase transitions features, and in particular the critical behavior found in 
ID lattice SP establishment, leads to modeling by a percolation process. Percolation theory 
belongs to a universality class that unifies the description of many physical models. As a 
direct consequence of the established properties of percolation, this framework allows us a 
new insight of SPs, without the approach implying exact solving of Maxwell's equations. 

The long- wavelength terahertz domain offers the advantage of the precise control of the 
geometry 16| . Therefore, the filling parameter of the subwavelength structures is accurately 
calibrated. Furthermore, the realization of ultra-thin (A/50) metal sheets removes possible 

n 

wave-guiding Fabry-Perot resonances [17|]. The structures are made of free-standing electro- 
formed ultra-thin nickel plates, with precision of design better than 1 /an. Terahertz spectra 
are recorded using; terahertz t im e- dom a,n spectroscopy N . Broadband linear* po.ar.ed 
subpicosecond single cycle pulses of terahertz radiation are generated and coherently de- 
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tected by illuminating photoconductive antennas with two synchronized femtosecond laser 
pulses. The sample is positioned on a 10 mm circular aperture, in the linearly polarized, 
frequency independent, 4.8 mm-waist Gaussian THz beam (1/e in amplitude). Numerical 
Fourier transform of the time-domain signals with and without the sample gives access to 
the transmission spectrum of the subwavelength structures. The dynamics of the surface 
plasmons is recorded over 240 ps, yielding a 3 GHz frequency precision after numerical 
Fourier transform, with 10 4 signal to noise ratio in 300 ms acquisition time. 

We studied the generation of SPs by analyzing the terahertz transmission spectra of 
series of free-standing subwavelength ID and 2D lattices (figures HJA.-C) with increasing 
quantity of metal, while the periodicity remains constant. A filling parameter p is defined 
as the ratio between the surface of the metal S m and the surface of the apertures S a in 
the lattice, by p = S m j S a , which enables to compare lattices of different shapes. SPs are 
built when increasing this filling parameter up to a precise level. Transmission spectra 
show typical resonant features, with minimum in ID lattice (figure HP), maximum in 2D 
lattice (figure [HE), and for both, a frequency shift Av = v — Uq of the first resonance 
is associated with the metal filling, where uq is the constant resonance frequency at low 
fillings. Indeed, in thin metal grids, resonant profiles due to diffraction over the periodic 
structure are observable near the transition and still allow defining a resonance frequency. 
Furthermore, we define the energy U as the total electromagnetic energy blocked by the 
plate, and can be directly accessed from the transmission of the terahertz beam through the 
system U = 1 — f \A t (t)\ 2 dtj f \A (t)\ 2 dt, where A t and A are the measured amplitudes 
of the electric fields with and without the sample, respectively. We then define a geometric 
capacity C linking the internal energy to the filling parameter, by C = dU / dp. We focused 
on the observation and evolution of these parameters, which provide fundamental insights 
on the nature of the generation of SPs. 

Series of two-dimensional arrays of squared and round subwavelength holes (36 lattices 
each) have been analyzed (figure dB-C). Incident linear polarization is sent along direction 
x. At low metal fillings, transmission spectra (figure [TJ5 black) show two large resonance 
features with a factor 2 between the resonance frequencies and no frequency dependent 
shift. The incident electromagnetic wave is diffracted through the subwavelength apertures 
(Bethe's theory [19]) and the transmitted waves interfere over the hole array. At high values 
of the filling parameter on the contrary, typical features of 2D lattice SPs can be observed 
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(figure [JJE color) with sharp maxima and resonance frequencies scaling as \/2. A resonance 
frequency shift with p is also observed, with precision better than 2 GHz. For tiny holes, first 
and second resonances agree well with theoretical resonance frequencies given by Bloch wave 
model at 0.5 and 0.71 THz, respectively |4(, for a lattice period of L = 600 fim. Increasing the 
filling parameter, experiments show clear discontinuity of the first resonance frequency and 
geometric capacity, at a specific transition filling parameter, both for squared (figure [2]A. and 
C) and round holes (figure [2)3 and D). This behavior is highly evocative of a first order phase 
transition. Transition filling parameter is the same for frequency shift and capacity, but 
different between squared and round holes, p = 1.03 ±0.01 and p = 1.37 ±0.01, respectively. 
It may also be related to a cut-off for the polariton propagation along the surface. Frequency 
shift is constant below the transition filling parameter and then monotonously increases after 
the discontinuity in the SP phase. It converges toward the value given for tiny holes, which 
is the same for squared and round holes [3J]. Frequency shift and capacity show different 
behavior whether the holes are squared or circular, with different transition filling parameter 
and discontinuity amplitude, confirming that SP is strongly influenced by the shape of the 
subwavelength holes. In addition, we have performed numerical simulations. We carried out 
a direct resolution of Maxwells equations through a full three-dimensional (3D) ab initio 
finite Element Method (FEM) analysis of the electric field propagating trough the array 
is) ] . Results of the frequency shift for square and round apertures are shown in figure [2], 
and demonstrate clear discontinuity at the same filing parameters found in the experiments. 

The one-dimensional subwavelength lattices are made of very narrow slits (figure [DA.). 
Series of ID subwavelength lattices of increasing filling has been realized with constant lattice 
period L =400 /iin. The terahertz beam is sent perpendicularly to the lattice surface and 
linear polarization has an angle 9 with respect to the transverse direction y. At low fillings, 
the subwavelength slits behave as standard polarizers. In particular, the spectra show a 
modulation with minimum resonances given by the diffraction over the metal lines (figure [Tp 
black), and the amplitude spectra scale proportionally to cos(6'), with no resonance shift. 
On the contrary, at high fillings, the spectra show evidence of typical ID SPs (figure [Tp red) 
with nearly zero minimum resonance of transmission (0.75 THz at 9 — 0°). Furthermore, 
this minimum shifts toward the low frequencies when incident polarization is rotated, due to 
the change of the projected apparent period. We measured the frequency shift for a set of 60 
lattices of different filling parameters. For each angle 9, the evolution of the frequency shift 
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with the filling is characteristic of a second order phase transition (figure [3]A). The frequency 
shift remains constant in the low filling range, up to a critical value p c = 1.583 ± 0.005, from 
which it continuously increases. This characteristic behavior is similar for all the incident 
polarizations. Evolution of the frequency shift versus the incident polarization angle shows 
sinusoidal variations (figure [3]B) as Au = Ausm(9 — 9o), where Av is the intrinsic frequency 
shift and 8q a possible small experimental angular misalignment. Precision on Au is better 
than 1 GHz at maximum deviation. Defining a normalized filling parameter r = {p — p c )/Pc, 
the evolution of the frequency shift in the SP domain near p c follows Au oc |r| , where (3 
is the critical coefficient associated with the frequency shift. Using Au, one obtains the 
critical coefficient (3 = 0.98 ± 0.04 (figure Geometric capacity is also analyzed from 
internal energy. Capacity shows diverging features near p c (figure H]) . Geometric capacity 
evolves as C oc |r| a , where a is the critical coefficient associated with capacity. For both 
sides of p c (see inset in figure H]), this critical coefficient has a value of a + = —0.99 ± 0.08 
and ct_ = —0.97 ± 0.08, above and below p c , respectively. 

The analysis of the experimental data is now discussed. The main first point is that 
the data can not be fully understood by existing models of SPs. These models intrinsically 
consider the existence of only one possible state for the surface wave. Waveguide analysis 
or Bloch wave analysis provide an approximation of the resonance frequency 2(]] and give 
good values only in the limit case where the apertures are very tiny. Further refinements 
consist in decomposition on a basis, such as Rayleigh expansion, but their linear features 
fundamentally prohibit discontinuity. This discrepancy is illustrated from 21[ in figure [2] 
(dashed lines). On the contrary, we have observed here a first order transition behavior 
in 2D lattice, associated with the simultaneous discontinuity of both frequency shift and 
geometric capacity at the same filling value. Moreover, in ID lattice, power law evolution 
of frequency shift and capacity is intimately linked to second order phase transition, and 
the critical coefficients are clearly in favor of a percolation process. Another strong hint is 
the polarization response below and above p c . There is a spontaneous symmetry breaking 
phenomenon at p c , since original isotropic response below p c is broken. Indeed, at low fillings, 
the shape of the spectra remains unchanged with the rotation of the incident polarization. 
On the contrary, the frequency resonance shifts with the polarization angle above p c (see 
figure [3]A.) , corresponding to a disappearing of the symmetry. Therefore, the generation of 
SPs can be described in a phase transition framework. The first phase (p < p c ) consists of 
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non resonant surface waves diffracted by the lattice in a way similar to Bethe description 
[3]. The second phase (p > p c ) is SPs characterized by a resonant interaction with the 
lattice. 

Complex systems can be described by statistical physics, in particular systems involving 
geometrical consideration, such as percolation in conductivity, polymerization or forest fire 
propagation. In the same way, in ferromagnetic/paramagnetic transitions, only the measure- 
ment of the linear magnetic susceptibility gives access to the second order phase transition. 
Phase transitions imply local non linear interactions. Even though the measured quantities 
are linear with respect to the electromagnetic wave going through the plate, the variation of 
the size of the aperture changes the relative contribution of far- and near-field interactions, 
which behave as 1/r and 1/r 3 . Thus, non linear features with respect to the geometry are 
the key of the understanding of our system. Furthermore, the screen transmission is not 
linear with the number of holes, since the transmitted signal is not the sum of the signal of 
a single hole. 

From macroscopic measurements, we demonstrated the phase transition origin of SP gen- 
eration. Universality class theory shows that the critical exponents for the phase transitions 
for broad classes of physical systems are the same. Therefore, a model describing a specific 
system can be extended to another system within its universality class. Then, experimental 
critical values of a and (3, associated with universality of percolation theory provide us a 
microscopic description of SP generation. The first order phase transition is also described 
by percolation and the difference between ID and 2D lattices originate from the change of 
degree of freedom in the system. Percolation transition describes microscopic "sites" linked 
together by "bonds". Consequences of percolation theory are that in SPs the sites are the 
(randomly moving) surface electrons, and the bonds are the coupling between the electrons 
and photons from the electromagnetic wave. Surface electrons are locally correlated by the 
interaction with the electromagnetic wave, defining a cluster of electrons as an ensemble 
of correlated surface electrons. A percolation transition then occurs when the first cluster 
reaches the size of the screen. At low filling factor, the size of the clusters is small, electrons 
only interact locally, and no SPs are present (cf figure [Tp-E). When p > p c , large groups of 
surface electrons interact together, forming clusters of macroscopic size and establishing a 
SP. Furthermore, the experimental critical parameters a ~ — 1 and (3 ~ 1 can be related to 
the percolation class of Cayley tree 22j |. which provides the critical evolution of the cluster 
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size, as x oc |p — Pep 7 , with 7 = 1. 

It has been demonstrated that SPs are a phase of surface electrons coupled to photons, 
and that SP building implies percolation phase transitions with respect to a filling parameter. 
In 2D lattice, a first order phase transition exhibits discontinuity of order parameter and 
geometric capacity, associated with a hole shape dependent transition filling parameter. In 
ID lattice, second order phase transition was analyzed and two critical coefficients measured, 
associated with Bethe lattice percolation model. Under critical conditions, correlation length 
diverges, with the consequence that information can be coherently transmitted through the 
plasmonic device on very large distance. Furthermore, ID SP lattice appears to be promising 
easy-to- handle experimental tools for studying second order phase transitions. 
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Figures 



FIG. 1: Geometry of the ID (A) and 2D (B-C) lattices. Terahertz transmission spectra. 

(D) ID lattice spectra, at low filling (p = 1.5, black line) and high filling (T=1.67, 6 = 0° solid 
red line, 6 = 80° dashed red line). (E) 2D lattice spectra of squared holes, at low filling {p = 0.78, 
black line) and high filling (p = 3.3 green line, p = 7.2 blue line, p = 9 red line). Arrows show 
theoretical Bloch wave theory frequencies. 



FIG. 2: First order phase transition in 2D lattices. Evolution of the frequency shift for 
squared holes (A) and round holes (B). Solid lines are 3D FEM numerical simulations of the 
evolution of the frequency shift. Evolution of geometric capacity for squared holes (C) and round 
holes (D). Transition filling parameters for squared and round holes are p = 1.03 ± 0.01 and 



p = 1.37 ± 0.01, respectively. Dashed lines are linear calculations from 



FIG. 3: Second order phase transition in ID lattice: resonance frequency shift. (A) 

Evolution of the frequency shift Av with filling parameter and corresponding slit width, for incident 
polarization angle =80° (black), 45° (red) and 20° (blue). Critical filling parameter is p c = 
1.583 ± 0.005. (B) Evolution of the frequency shift with incident polarization at p = 1.67 for 
experimental data (solid circles), and sinusoidal fit (solid line). (C) Evolution of the frequency 
shift versus normalized filling parameter in logarithm scale showing the critical behaviour for data 
(circles) and critical evolution fit with (5 = 0.98 ± 0.04 (solid line). 
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FIG. 4: Second order phase transition in ID lattice: geometric capacity. Evolution of 
capacity versus filling parameter for experimental data (circles) and critical divergence fit (solid 
lines). Inset shows critical evolution at both sides of the critical filling parameter, with a + = 
—0.99 ± 0.05 and ct_ = —0.97 ± 0.05, above and below p c , respectively. 
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Masson et al. Figured] 
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Masson et al. Figure [2] 
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